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I. INTRODUCTION 

Recently, there has been a growing interest in investi- 
gation of gyraton spacetimes. They represent the grav- 
itational field of a localized source with an intrinsic ro- 
tation, moving at the speed of light. Such an idealized 
ultrarelativistic source, which can be modeled as a pulse 
of a spinning radiation beam, is accompanied by a sand- 
wich or impulsive gravitational wave. 

In fact, gravitational fields generated by (nonrotating) 
light pulses and beams were already studied by Tolman 
1] in 1934, who obtained the corresponding solution in 
the linear approximation of the Einstein theory. Ex- 
act solutions of the Einstein-Maxwell equations for such 
"pencils of light" were found and analyzed by Peres 
and Bonnor [1, 0, IE] • These solutions belong to a general 
family of pp -waves [f| 0] • 

In the impulsive limit (i.e. for an infinitely small cross- 
section of the beam, and for the delta-type distribution 
of the light-pulse in time), the simplest of these solu- 
tions reduces to the well-known Aichelburg-Sexl metric 
Q which describes the field of a point-like null particle. 
It can be obtained by boosting the Schwarzschild metric 
to the speed of light, with the mass tending to zero so 
that the total energy is kept finite. More general impul- 
sive waves were subsequently obtained by boosting other 
black hole spacetimes with rotation, charge and a cos- 
mological constant 0, M EL HI HE 0, m (for recent 
reviews, see [iH ITm]). 

Gyraton solutions are special sandwich or impulsive 
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waves of the Kundt class (which generalize the pp -waves) 
such that the corresponding beam of radiation carries not 
only energy but also an additional angular momentum. 
Such spacetimes were first considered by Bonnor in [l8|, 
who studied the gravitational field created by a spinning 
null fluid. He called the corresponding particle made out 
of this continuum a "spinning nullicon" . In some cases, 
this may be interpreted as a massless neutrino field [lj|. 

In the exterior vacuum region outside the source, these 
solutions are locally isometric to standard pp -waves. 
The interior region contains a nonexpanding null matter 
which has an intrinsic spin. In general, these solutions 
are obtained by keeping the nondiagonal terms g U i in the 
Brinkmann form [20( of the pp-w&ve solution, where u is 
the null coordinate and x 1 are orthogonal spatial coordi- 
nates. The corresponding energy-momentum tensor thus 
also contains an extra nondiagonal term T u j = fa. In four 
dimensions, the terms g U i can be set to zero locally, using 
a suitable gauge transformation. However, they can not 
be globally removed because the gauge invariant contour 
integral § g u i(u, x 3 ) dx l around the position of the gyra- 
ton is proportional to the nonzero angular momentum 
density ji, which is nonvanishing. 

Similar gyratons in a higher dimensional flat space 
were investigated (in the linear approximation) by Frolov 
and Fursaev [2l| . Such gyratons represent a pulse of cir- 
cularly polarized radiation or a modulated beam of ul- 
trarelativistic particles with spin or other sources, which 
have finite energy E and finite total angular momen- 
tum J. The gyraton itself is characterized by two ar- 
bitrary profile functions of u which determine the energy 
density and angular momentum. The authors investi- 
gated the limit in which the source becomes infinitesi- 
mally small (with a negligible radius of the cross-section) 
and the profile functions are independent. They also 
studied the geodesic motion of test particles in the field of 
gyraton and demonstrated that, when the gyraton passes 
through the center of the ring of test particles, the par- 
ticles start to rotate. In fact, the gyraton's angular mo- 
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mentum effectively creates a force which is similar to the 
usual centrifugal repulsive force, while the gyraton en- 
ergy produces the attractive "Newtonian" force. 

Frolov, Israel, and Zelnikov [12] further investigated 
the exact gyraton solutions propagating in an asymp- 
totically flat D-dimensional spacetime and proved that 
the Einstein's equations for gyratons reduce to a set of 
linear equations in the Euclidean (D — 2)-dimensional 
space. They also showed that the gyraton metrics be- 
long to a class of vanishing scalar curvature invariants 
(VSI) spacetimes for which all polynomial scalar invari- 
ants, constructed from the curvature and its covariant 
derivatives, vanish identically [Hj]. (For the discussion of 
spacetimes with nonvanishing but nonpolynomial scalar 
invariants of curvature, see [2J].) Subsequently, charged 
gyratons in Minkowski space in any dimension were pre- 
sented in df|. 

In [111, Frolov and Zelnikov took a cosmological con- 
stant into account, and exact solution for gyratons in the 
asymptotically anti-de Sitter spacetime were presented. 
Namely, they obtained Siklos gyratons which generalize 
the Siklos family of nonexpanding waves [27j (investi- 
gated further in [H]). 

In this case, all polynomial scalar invariants are in- 
dependent of the arbitrary metric functions which char- 
acterize the gyraton and have the same values as the 
corresponding invariants of pure anti-de Sitter back- 
ground. The AdS gyratons [lH thus belong to the 
class of spacetimes with constant scalar invariants (CSI) 
[29l l30l l3ll l32l . [33j . In string theory it has been demon- 
strated that generalized pp-wave spacetimes do not get 
any quantum and a' corrections and hence arc perturba- 
tively exact. One may expect a similar property to be 
valid also for the gyratons, but more careful analysis is 
required since even if all of the local counterterms in the 
effective action are trivial constants for CSI spacetimes, 
their metric variations can be nontrivial functions. Still, 
one can try to generalize the property of relatively simple 
quantum corrections to the case of semiuniversal metrics 
[3fj ] when the Ricci tensor has a block-diagonal structure. 

Let us also mention that string gyratons in supergrav- 
ity were recently found in [34[. Supersymmetric gyraton 
solutions were also obtained for a minimal gauged theory 
in five dimensions in [3o| , where the configuration repre- 
sents a generalization of the Siklos waves with a nonzero 
angular momentum in anti-de Sitter space. 

The gravitational field generated by gyratons may be 
interesting for studies of production of mini black holes in 
colliders (such as the LHC) or in cosmic ray experiments. 
The problem of mini black hole formation in high energy 
particle collisions is an important issue of TeV gravity. 
The theory of such collisions, developed in J36l [371 1381 
[39l . lioj , was applied to gyraton models in [4lj] . 

The purpose of our contribution is to further extend 
the family of gyratonic solutions, which are only known 
in Minkowski or anti-de Sitter background spaces. In 
particular, we present a new large class of gyratons of al- 
gebraic type II, propagating in less trivial universes which 



are a direct product of two 2-spaces of constant curva- 
ture. This family of vacuum and electrovacuum back- 
ground spacetimes contains the Nariai [H|], anti-Nariai, 
and Plebahski-Hacyan universes 14311 of type D, or con- 
formally flat Bertotti-Robinson [4J, [45| and Minkowski 
spaces. These direct-product spacetimes with six isome- 
tries (see @, 0] for more details) recently attracted new 
interest because they can be recovered as specific extreme 
limits of various black hole spacetimes in four or more di- 
mensions Hi E3, m, m • 

Impulsive gravitational and pure radiation waves in the 
(anti-)Nariai, Bertotti-Robinson, and Plebahski-Hacyan 
universes were presented and analyzed by Ortaggio and 
Podolsky [H3, HH. They showed, and subsequently an- 
alyzed in more detail in that these solutions are 
straightforward impulsive limits of a more general class 
of Kundt spacetimes of type II with an arbitrary profile 
function, which can be interpreted as gravitational waves 
propagating on specific type D or O backgrounds, includ- 
ing those which are a direct product of two 2-spaces. In 
fact, the gyraton spacetimes investigated in this paper 
are generalizations of such Kundt waves when their ul- 
trarelativistic source is made of a "spinning matter" . 

The paper is organized as follows. In Section [Til we 
present the ansatz for the metric and fields. After a short 
review of the transverse space geometry we derive field 
equations and simplify them introducing the potentials. 
Next, we discuss the gauge freedom and suitable gauge 
fixings. The overview of the gyraton solutions is summa- 
rized in Section. Ill Fl 

In Section IIII1 we give a survey of important special 
subclasses of our gyraton solution. They include direct- 
product spacetimes, all type D vacuum backgrounds, and 
general Kundt waves on these backgrounds. In Section 
IIV1 we concentrate on the interpretation and description 
of the gyratons. We discuss geometric properties of the 
principal null congruence, the Newman-Penrose (NP) 
quantities with respect to natural tetrads, and proper- 
ties of the electromagnetic field. 

The final Section V describes the Green functions re- 
quired to solve the field equations. The main results of 
the paper are summarized in concluding Section [VII 

Some technical results needed to derive the field equa- 
tions and NP quantities are left to Appendices [5] and [B] 
In Appendix [Ul we derive all electro- vacuum solutions of 
type D, and Appendix [D] discusses further details con- 
cerning the Green functions. 



II. GYRATONS ON DIRECT PRODUCT 
SPACETIMES 

A. The ansatz for the metric and matter 

The aim of this paper is to derive and analyze the 
family of gyraton solutions describing a gyratonic matter 
which propagates, together with a related gravitational 
wave, through a direct-product spacetime filled with a 
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"uniform" electromagnetic field. 

We assume that such spacetimes belong to the Kundt 
class. It is characterized by a geometrical property that 
it admits a nonexpanding, nontwisting, and shear-free 
null congruence [6j, l7J. This congruence represents the 
null direction of propagation of the gyraton and of the 
accompanying gravitational wave. 

In terms of canonical (real) coordinates {r, u, x,y}, 
such a metric reads 

ds 2 =7^ (da; 2 + dy 2 ) - 2dudr - 2H du 2 

(2.1) 

+ 2a x dx du + 2a y dy du , 

where H(r,u, x,y) can depend on all coordinates, but 
the functions a x (u, x, y), a y {u,x,y), and P{u,x,y) are 
r-independent. The restriction d r P = follows from our 
assumption of vanishing expansion of the Kundt geome- 
try, while the condition d r ai = 0, where i — x,y, is nec- 
essary here to obtain a gyraton which propagates on a 
direct-product spacetime background. In fact, this con- 
dition is a consequence of the Maxwell equations in the 
case when the electromagnetic field is present. In the ab- 
sence of the electromagnetic field, the vacuum Einstein 
equations admit that functions dj can be linear in r. 
However, geometrical properties of such solutions are 
substantially different from those of the direct-product 
spacetimes. Therefore, in the following we will always 
assume that 

d r a x = , d r a y = . (2.2) 

This assumption thus implies that such solutions belong 
to the special subclass of Kundt solutions (see section 
HVBI for more details). 

The metric should satisfy the Einstein equations with 
a stress-energy tensor generated by the electromagnetic 
field and the gyraton: 

G MV + Ag uv = x(T™ + 1*7) . (2.3) 

Here, A and x = 8ttG are the cosmological and gravita- 
tional constants, respectively. 

The spacetime can be filled with the background elec- 
tromagnetic field, which is modified by a gravitational 
influence of the gyraton. We assume 

F = E dr A du + B dx A dy 

+ a x du A dx + a y du A dy , 

where E and B are constants, so that the corresponding 
stress-energy tensor T™ has the form (|A3|) . This ansatz 
for the Maxwell tensor has been inspired by the electro- 
magnetic field known in the Bertotti-Robinson [44], [4f| 
and Plebanski-Hacyan spacetimes [43| , to which we have 
added new terms proportional to functions <Ji(r,u,x,y). 
In fact, terms with such a structure are generated if we 
demand a gauge symmetry of the electromagnetic field 
under gauge transformation discussed in III El 



Finally, we must characterize the gyratonic matter by 
specifying the structure of its stress-energy tensor. It is 
a generalization of a standard null fluid such that we ad- 
ditionally allow terms corresponding to 'internal spatial 
rotation' of the gyraton source, 1 

x T syl = j u du 2 + 2j x du dx + 2j y du dy . (2.5) 

We admit a general coordinate dependence of the source 
functions u, x, y) and ji(r, u, x, y). However, it will 
be shown below that the field equations enforce a rather 
trivial r-dependence of these functions. Let us note that 
previous papers on gyratons, namely [22j], [HI, (25[, as- 
sumed that the gyraton source is r-independent. 

The gyraton source is thus described only on a phe- 
nomenological level, by its stress-energy tensor (|2.5[) . We 
do not discuss a possible internal structure of the gyra- 
tonic matter, and we do not specify its own field equa- 
tions. The gyraton stress-energy tensor is assumed to be 
given, and our aim here is to determine its influence on 
the metric and the electromagnetic field. However, we 
have to consider that the gyraton stress-energy tensor is 
locally conserved. It means that the functions j u and % 
must satisfy the constraint given by 

T~ r = . (2.6) 

Of course, if we had considered a specific internal struc- 
ture of the gyratonic matter, the local energy-momentum 
conservation would have been a consequence of field equa- 
tions for the gyraton. Without that, we have to require 
(123)) explicitly. 

To summarize, the fields are characterized by functions 
P, H, di, and <Ji which must be determined by the field 
equations, provided the gyraton sources j u and ji and the 
constants E and B of the background electromagnetic 
field are prescribed. 

As we will discuss in IIII Al pure background solu- 
tions are obtained when both gyratons and the gravi- 
tational waves are absent, namely for T^fJ = 0, a* = 
and Ft oc r 2 . For the Minkowski and (anti-)Nariai back- 
grounds, also vanishes, while it is nonzero for the 
Bertotti-Robinson and Plebanski-Hacyan spacetimes. 

Finally, for later convenience, we introduce a con- 
stant p, given by the parameters E and B of the elec- 
tromagnetic field, 

p= ^(E 2 + B 2 ), (2.7) 

(with k and e being gravitational and electromagnetic 
interaction constants, 2 respectively) and the constants 
A + and A_, defined as 

A± = A±p. (2.8) 



1 In all tensorial expressions for the metric and other symmetric 
tensors, we understand by, for example, dudx the symmetric 
tensor product | (du ® dx + dx (g> du) . 

2 There are two natural choices of geometrical units: the Gaussian 
with x = 87T, and e Q = 1/An, and Si-like with h = e D = !• 
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B. Geometry of the transverse space 

The geometrical structure of the Kundt metric (|2.ip 
identifies the null geodesic congruence generated by d r 
and parametrized by an affine time r, the family of null 
hypersurfaces u — constant, and two-dimensional trans- 
verse spaces r,u = constant. It will be convenient to re- 
strict various equations to these transverse spaces. For 
example, <2j and Oi can be understood as components 
of u-dependent 1-forms on these two-dimensional spaces. 
Therefore, we now briefly review some formulas and def- 
initions valid in such two-dimensional transverse geome- 
try. 

The transverse space is covered by two spatial coordi- 
nates x 1 , and we use the Latin indices i, j, . . . to label the 
corresponding tensor components. The restriction of the 
metric (|2.ip to the transverse space is 

ds 2 ± = g ±ij dxW = (dx 2 + dy 2 ) . (2.9) 

Here we made a useful choice of coordinates x % = {x, y} 
in which ds 2 ^ has a conformally flat form. 3 

The transverse curvature is fully characterized by the 
scalar curvature R±, which in terms of conformally flat 
coordinates reads (cf. the definition (|2.18p below) 

^R ± = A lOgP = />(/>., + Pyy) - {P% + P 2 y ) . (2.10) 

Inspecting the ru component of the Einstein equa- 
tions (|2.3p . we find that the transverse scalar curvature 
has to be constant, 

\Rx = AlogP = A+, (2.11) 

cf. the first lines in equations (|A3[) and (|A4p . together 
with (|2.8p . The transverse spaces are thus the constant 
curvature 2-spaces, all with the same curvature. Thanks 
to this property we can further simplify the choice of the 
transverse coordinates {x, y} in such a way that the con- 
formal factor P~ 2 in (|2.9p is u- independent. Therefore, 
in the following we may assume 

d r P = , d u P = . (2.12) 

Moreover, using a freedom in the choice of the trans- 
verse coordinates, we can also put the conformal factor 
P to a canonical form. There are two standard choices 
solving (|2.1ip . namely, 

P=l + U + (x 2 + y 2 ), (2.13) 



3 The conformally flat coordinates are not essential, but they sim- 
plify some expressions. In a two-dimensional space, a choice of 
such coordinates is always possible. 



and, for a negative A + , 

P = V-A+ x . (2.14) 

However, in the following, we do not need a particular 
form of P. It must just satisfy Eq. (|2.11l) . 

With the transverse metric l|2.9p we may associate the 
Levi-Civita tensor (with e xy = P~ 2 ) and the covariant 
derivative denoted by a colon (e.g., di-.j). We raise and 
lower the Latin indices using g ±i j , and we use a shorthand 
a 2 = a l ai — P 2 (a 2 + a 2 ) for a square of the norm of a 
1-form a,. In two dimensions, the Hodge duals of 0-,l- 
and 2-forms ip, a,-, and read 

(*(p)ij = (peij , (*a)i = aje 3 , , */ = ^fijf? 3 = P 2 fx V ■ 

(2.15) 

For convenience, we also introduce an explicit notation 
for two-dimensional divergence and rotation of a trans- 
verse 1-form di, 

diva = a{ % = P 2 {a x<x + a y>y ) , ( 2 -16) 
rota = *da = = P 2 (a VyX — a x . y ) , (2-17) 

and for the Laplace operator of a function ijj, 

A^ = ^ 1 =P 2 (^, xx +^,yy) . (2.18) 

Note that the divergence and rotation are related as 
div a = rot *a. 

Finally, we will generally assume that the transverse 
space is topologically simple in the sense that the space 
of harmonics is trivial. However, sometimes it will be 
physically relevant to consider also nontrivial solutions 
of the Laplace equation if we relax the boundary and 
asymptotical conditions in the noncompact case. For ex- 
ample, a solution of the Laplace equation around a local- 
ized source satisfies the homogeneous Laplace equation 
on the space with the source removed. Such a space is, 
however, noncompact and the solution is not vanishing 
on the boundary. 

C. The field equations 

After specifying the ansatz for our fields and reviewing 
the transverse geometry we can now derive the equations 
for the gyraton. We have to consider the Einstein equa- 
tions (|2.3p together with the Maxwell equations and the 
condition (|2.6p for the gyraton source. 

We start with the cyclic Maxwell equation. Assuming 
(gill) and (j2"T2"|) it reads 

= dF = (d r (Ti)dr AdwA da; 1 — rot adu Ae , (2.19) 

where e = P~ 2 dx A dy. We immediately infer that the 
1-form <Xj is r-independent, d r Ui = 0, and rotation-free, 

rotcr = 0. (2.20) 
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The second Maxwell equation F av ' v = has only 
the u component nonvanishing, 4 which gives 

div a — E div a + B rot a = . 



(2.21) 



We call (|2.20p and (|2.2ip the potential equations since 
they guarantee the existence of potentials which will be 
discussed in detail in Section TlIDI For this reason, it is 
useful to note that these equations imply the conditions 



div [E {a-E a) + B *(a-E a)] = , 
rot [E (cr-B *a) + B *(a~B *o)] = , 



(2.22) 



cf. equations (|2.39[) below. 

The Einstein equations can be derived from the Ein- 
stein tensor and the electromagnetic stress-energy ten- 
sor, which are given in Appendix [X] We have already 
discussed the ru component which leads to the condi- 
tion (|2.1ip . The transverse diagonal components xx and 
yy give 



dlH 



-A_ 



(2.23) 



We thus obtain the explicit r dependence of the metric 
function H as 



gr 



(2.24) 



where we have introduced r-independent functions 
g{u,x^) and h^u^x^). 

Finally, the remaining nontrivial components of the 
Einstein equations are those involving the gyraton source 
(|2.5p . The ui components give an equation related to ji , 
which we call the first source equation, 



ji = - fij' 3 +9,i - A- </, 



+ xe [E (o-i-Ea,i) + B (o-j-Eaj) e J i] 
I ' (2-25) 

+ xe [E (<Ti-B ajei i) + B (o-j-B a k e k j)e j i] , 

where we have introduced the external derivative fij of 
the 1-form a, as 



*3,i 



■ rot a) 



(2.26) 



For convenience, we have written the equation (|2.25p in 
two equivalent forms. In the square brackets, they explic- 
itly contain the terms which were already encountered in 
the equation (|2.22p . We can thus easily split the first 
source equation into divergence and rotation parts: 

-div j = -A g + A_ diva , (2.27) 
-rotj = -A6 + A+6, (2.28) 



where the function 6(m,x : ') is the Hodge dual of /y, 

b = *f = rota. (2.29) 



Equations (|2.27p and (|2.28p carry essentially the same 
information as the original source equation (|2.25|) . 5 

Next, we examine the condition (|2.6p for the gyraton 
source. It gives 

- {d r ji) dx l + (-drju + div j + a'drji) du = 0, (2.30) 

so that the source functions ji must be r-independent 
and j u has to have the structure 



j u = rdiv j + i 



(2.31) 



The gyraton source (|2.5p is thus fully determined by three 
r-independent functions l{u, x j ) and ji(u, x^). 

Finally, from the MM-component of the Einstein equa- 
tion we obtain 



j u = (a g — A_div a) r 

+ Ah + ^b 2 - A_a 2 + 20*3,., 

+ d u (div a) + g div a 
— xe Q (a — Ea) 2 . 



(2.32) 



Comparing the coefficient in front of r with (|2.2T[) . we 
find that it consistently reproduces the structure (|2.3ip . 
The nontrivial r-independent part of (|2.32[) gives the 
second source equation which can be understood as the 
equation for the metric function h, 



A h = l - -b 2 + A-a 2 - 2a?g i 

2 y ' 1 (2.33) 

+ yce a {a — Ea) 2 — (9„(diva) — gdiva . 



D. Potentials 

We have thus found that the Maxwell and Ein- 
stein equations reduce to two potential equations (|2.20[> . 
(|X2"I|) . and two source equations 1[2~25]) . (|2T52]I . These 
equations can further be considerably simplified by in- 
troducing potentials for the 1-forms Cj, Hj and for the 
source ji . 

Indeed, the first potential equation (|2.20p gives imme- 
diately that (Ji has a potential (p(u,x J ) such that 



(2.34) 



Using the Hodge decomposition we can express the 1- 
form ai using two scalar potentials k(u,x j ) and X(u,x J ): 



= Ki + et 3 A. 



(2.35) 



4 Here we used that and Cj are r-independent. If the condition 

(I2.2li was not assumed before, it would follow from the transverse 5 They are equivalent to l|2.25|l if we ignore the possibility of har- 
components of this Maxwell equation. monic 1-forms which can exist in topologically nontrivial spaces. 
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These potentials control the divergence and the rotation 
of <Zj via 

diva = AK, rota = — AA. (2.36) 

Equation (|2.21[) imposes a constraint among these 
three potentials tp, K, and A: 

A tp = A (Ek + BX) . (2.37) 

If the transverse space is compact (or if it is noncom- 
pact but sufficiently strong asymptotic conditions are im- 
posed) the solution of the Laplace equation is trivial and 
we immediately obtain 

tp = En + BX . (2.38) 

By using this constraint, it is possible to show that the 
potentials k and A solve the conditions (|2.22[) 

E {cn-Edi) + B fa-Eaj) e 3 \ = (E 2 + B 2 ) A. ,< ? , , 

E {(T.-B^ai) + B (aj-B*aj) e j i = [E 2 + B 2 ) K ti . 

(2.39) 

In terms of the potentials, the first source equation 
(I2.25P can be written as 

Ji = 2 fa '' j + A + A j eJ i - A - K ,i + 9,i ■ ( 2 - 4 0) 

Its rotation part is the equation (|2.28p for 6, the solution 
of which can be used as a source for the equation for the 
potential A, 

A\ = -b. (2.41) 

The divergence part of (|2.40|) can be written as a relation 
between the functions g, n and source divj, 

A(g- A_k) = divj . (2.42) 

The problem further simplifies if we introduce scalar 
potentials p(u, x 3 ) and q{u, x 3 ) for the gyraton source ji, 

ji = P,i + «/ qj , (2.43) 

so that 

rot j = — Ag, divj = Ap. (2.44) 

Substituting this to the field equation (|2.40|) and split- 
ting it into the gradient part and the rotation part (i.e., 
using the Hodge decomposition), we obtain 

g-K_K=p, (2.45) 

and 

^A A + A + A = -q . (2.46) 

Let us note that all of the potentials are defined up to an 
additive constant (which, however, can be u-dependent). 
In the derivation of (|2.45|) and (|2.46[) . we have absorbed 



the integration constants into this nonuniqueness of po- 
tentials. In view of (|2.41j) , function b is then given by 

b = 2(A + X + q). (2.47) 

We have thus reduced the field equations to simple al- 
gebraical relations (|2.38|) . (|2.45p between the potentials, 
to the Hclmholtz-Poisson equation (|2.46[) for A, and the 
Poisson equation (|2.33|) for h. The last one can be also 
rewritten using the potentials as 

Ah = c + qAX-pAK- 2a i p ii , (2.48) 

with h closely related to h: 

h = h + d u n+ ^A_k 2 - ^A+A 2 . (2.49) 

E. Gauge transformation and the field equations in 
suitable gauges 

Shift of the r coordinate 

To find the gyraton solution explicitly, we need to de- 
termine the functions h, g, a,i, and cr,, provided the gy- 
raton sources ji and i are prescribed. In terms of the 
potentials k, A, and tp, replacing the transverse 1-forms 
a, and cr,, we have obtained equation (|2.46[) for A, (|2.38p 
for tp, and (|2.33p for h. However, we have only one equa- 
tion (|2.42p for k and g. 

This deficiency of equations corresponds to the 
fact that our ansatz (j2~Tj) . I|2.4p . (|2.5p admits a 
gauge freedom. Indeed, the coordinate transformation 
f — > r — r — ip{u, x 3 ), accompanied by the following re- 
definition of the metric functions and fields: 

r = r — ip , 

g = g — A-ip , h = h- \^A^ 2 + g -0 + d u ip , 

ai = ai — ip t i, ai = ai— Eip ti , (2.50) 

k = k — ip , A = A, tp = (p — Eip , 

3i=3i , l = l + tpdivj , 

leaves the metric, the Maxwell tensor, and the gyraton 
stress-energy tensor in the same form. Consequently, all 
of the field equations remain the same. Such a transfor- 
mation is a pure gauge transformation and we can use it 
to simplify the solution of the equations. 

This gauge transformation has a geometrical meaning 
of shifting the origin of the affine parameter r of the null 
congruence d r . 

Inspecting this gauge transformation, we find that the 
combination g — A_k is gauge invariant. This combina- 
tion enters the field equation (|2.42p . and only this com- 
bination is thus invariantly determined by the sources, 
namely, it is equal to p, cf. (|2.45[) . The particular split- 
ting into g and k parts is just a question of the gauge 
choice. 
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Indeed, it follows from (|2.50|) that it is possible to mod- 
ify one of the functions g, k, or ip to an arbitrary value 
or even to cancel it out from all the equations. More- 
over, the freedom to choose one of these functions covers 
the gauge freedom fully. Therefore, we use them to con- 
trol the gauge freedom: we may fix the gauge by setting 
g, k, or ip to be an arbitrarily chosen function. Any of 
these gauge conditions leads to the same family of solu- 
tions, only with a different parametrization of the gauge 
freedom. 



Gauge fixing of g 



implies simple relations for g: 

Ag = divj , 
9=P> 

and between the 1-forms et^ and <7i (or their potentials): 



(2.57) 
(2.58) 



p = BX, 
Ui = B *at . 

The equation for h now takes the form 



(2.59) 
(2.60) 



Ah = t-2(A + X+q) 2 +A + X, i X, j g i: >+2X,iP, j e ij . (2.61) 



Let us start with the gauge condition that the func- 
tion g is an arbitrary function. Then the equation (|2.42p 
should be understood as the Poisson equation for n. In 
terms of the source potentials we have even the explicit 
solution given by (|2.45|) . The potential A is determined 
by the equation (|2.46[) . and tp by (|2.38|) . Finally, sub- 
stituting these results into the second source equation 
(|2.33[) we obtain the Poisson equation for h. 

In fact, the gauge fixing of g can be used to eliminate 
the metric function g completely Setting 



5 = 0, 

the field equation for k reduces to 

1 



A K : 



A. 



divj , 



with the solution 



- = -—_P 



(2.51) 



(2.52) 



(2.53) 



The equation for h reads 

Ah = i - 2(A + A + q) 2 +A + A., A ,.(/'•*' + 2 X ti pj e lj 



+ j-(p,iP,j9 ij +d u Ap) , 



(2.54) 



and the metric function H has only trivial quadratic de- 
pendence on r, 

H(r,u,x j ) = -^A_r 2 + h(u,x J ) . (2.55) 



Gauge fixing of « 

Alternatively, we can fix the potential k, which is 
equivalent to the prescription of a value of div a. Equa- 
tion (|2.42p is then the Poisson equation for g, otherwise 
the solution of the field equations proceed in the same 
way as above. 

The special choice 



Gauge fixing of tp 

The last natural gauge condition is a fixing of the po- 
tential (p. In this case, one first finds the potential A 
by solving the equation (|2.46[) . The relation (|2.38[) then 
gives the potential k. Plugging this into (|2.42[) , the equa- 
tion for g is obtained: 

Ag = divj + —rA<p-^-A\, (2.62) 



A_ A_J5 , 
9 = P+ -pr V? — A • 



E T E 

The 1-form a* can be written in terms of a; and A as 



(2.63) 



E a,i = <Ji — B X i — E Xj e 3 i 
Particularly, for 

tp = , i.e. o"j = , 

we obtain 

«=-§ A, 

g =p - A-— A , 
Ecu = -BXj - EXj^i . 



(2.64) 
(2.65) 

(2.66) 

(2.67) 
(2.68) 



The choice <p = simplifies the Maxwell tensor (12. 4p to 



F = E dr A du + B —dx A dy 
pi 



(2.69) 







i.e. div a — 



(2.56) 



Reparametrization of the u coordinate 

After the above discussion of the gauge freedom cor- 
responding to the transverse-dependent shift of the r co- 
ordinate, we should also mention the remaining gauge 
freedom. The metric (|2.1j) . the electromagnetic field 
(|2.4p , and the gyraton stress-energy tensor (|2.5[) keep the 
same form under a general reparametrization of the u 



8 



coordinate u — > u — f(u), accompanied by the rescaling 
f — > r = f/f'(u) of the r coordinate. The metric func- 
tions and matter fields must be redefined as 



= /(i 
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ft; 



/'(«) 



h = 



f'{u) 2 



J I 



/'(«) ' 

ju 
f'(u) 2 



A = 



o% = 



A 



I 



(2.70) 
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/'(«) ' 
Ji 



It is worth to emphasize here that this reparametrization 
is independent of the transverse spatial coordinates. This 
gauge transformation is thus 'global' from the point of 
view of the transverse space, and it does not influence 
the field equations (which we formulated as differential 
equations on the transverse space) in any significant way. 



The transverse 1- forms Of, <Tj, and ji can be written in 
terms of the scalar potentials «, A, <p, p, and q as 



°i = <P,i j 
d{ = K i -\- 6jp A j , 

ji = P.i + £t 3 Q,j ■ 



(2.77) 
(2.78) 
(2.79) 



These potentials are unique up to (unphysical) constants 
on the transverse space, which can always be gauged 
away. 

Finally, the functions h, g, k, A, and tp must satisfy the 
linear field equations 



p = Ek + BX , 

g- A_k = p , 



(2.80) 
(2.81) 



A + A + A = —q , 



(2.82) 

Ah = t + qAA-pAK - 2a i p >i , (2.83) 



where 



1 9 1 9 

h = h + d u n + - A-k z - 2 A + A 



(2.84) 



F. Summary of the gyraton solutions 



III. IMPORTANT SPECIAL SUBCLASSES 



Let us now summarize the main equations of the gy- 
ratons on direct-product backgrounds. These are space- 
times with the metric of the form 

ds 2 = i (dx 2 +dy 2 ) - 2 du dr - 2H du 2 + 2a t dx l du , 
pz 

(2.71) 

filled with the electromagnetic field 



F = E dr A du + B -^dx A dy + du A a l dx 1 , (2.72) 



and the gyratonic matter 

xT ! " = j u du 2 + 2j t dx l du . (2.73) 
The metric function H(r,u, x l ) is quadratic in r, 
1 



H = A_ r l + qr + h 

2 



(2.74) 



the gyraton energy density j u (r, u, x 1 ) can be at most 
linear in r, 



j u = r div j + i, 



(2.75) 



and the functions g{u,x l ), h(u,x l ), dj(u,x l ), <jj{u,x l ), 
jj(u,x l ), and l(u, x r ) are r-independent. The function 
P{x l ) is r and u-independent and it satisfies the equation 



A log P = A + 



(2.76) 



It can be solved by P of the form (|2.13[) or, for A + < 
by (EH. 



The large family of solutions of Einstein-Maxwell 
equations discussed above belongs to the Kundt class 
(12. ip of nonexpanding, shear-free and twist-free space- 
times [(| I531 . |54] , namely to its subclass characterized by 
the condition ()2.2|1 . As we have seen in Section Hi A[ and 
will be discuss more in Section ITVl the gyratonic matter 
(|2.5p is the "rotating" generalization of a null fluid. As 
special cases, this family of solutions contains some pre- 
viously known spacetimes from the Kundt family which 
correspond to electro- vacuum or pure (null) radiation. In 
this section we will shortly discuss such important sub- 
cases. 



A. Direct-product background spacetimes 

It is natural to start with the simplest case of highly 
symmetric spacetimes. Considering the metric function 
H of the form (|2.24[) . setting = 0, and g = = h, and 
choosing the expression (|2.13[) for P, the metric (|2.ip 
reduces to 



ds 2 = 



dx 2 + dy 2 



[1 + U+ (:r 2 + y 2 



2dudr + A_r 2 du 2 . (3.1) 



It describes backgrounds on which the gyratons 
propagate. By performing the transformation 



= v(l - iA_ uvy 1 with u=(t-z)/V2, 
v = (t + z)/\/2, the metric becomes 

dx 2 + dy 2 dz 2 



and 



ds 1 



dt 2 



[l + iA+^+y 2 )] 2 [l + |A_(z 2 -< 2 



(3.2) 
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Table I: Some of possible background spacetimes which are 
the direct product of two 2-spaces of constant curvature. Here 
A is a cosmological constant, and p is a constant energy den- 
sity of the electromagnetic field. 



A+ 


A_ 


geometry 


spacetime 


A 


p 








E 2 x M 2 


Minkowski 


= 


= 


A 


A 


S 2 x dS 2 


Nariai 


> 


= 


A 


A 


H 2 x AdS 2 


anti-Nariai 


< 


= 


P 


~P 


S 2 x AdS 2 


Bertotti-Robinson 


= 


> 


2A 





S 2 x M 2 


Plebanski-Hacyan 


> 


= A 





2A 


E 2 x AdS 2 


Plebanski-Hacyan 


< 


= |A 



The background spacetimes thus have geometry of a di- 
rect product of two 2-spaces of constant curvature A + 
and A_, respectively. The first is the space spanned 
by two spatial coordinates so that it is flat Euclidean 
space E 2 , 2-sphere S 2 , or 2-hyperboloid H 2 , according 
to the sign of the constant A+. The second is the (1+1)- 
dimensional spacetime spanned by a timelike coordinate 
and one spatial coordinate. According to the sign of the 
constant A_, it is Minkowski 2-space M2, de Sitter space 
dS-2, or anti-de Sitter space AdS%. 

Therefore, there are nine theoretically possible distinct 
subclasses given by the choice of A + and A_, but only 
six of them are physically relevant because the energy 
density p must be non-negative, which eliminates three 
cases. The most important of such background space- 
times are summarized in Table |H In addition, there are 
more general Bertotti-Robinson direct product space- 
times for which the constants A + and A_ are independent 
and nontrivial; i.e., the cosmological constant A and the 
energy density p > of the electromagnetic field can be 
chosen arbitrarily. 

In a natural null tetrad, the only nonvanishing NP 
Weyl and curvature scalars are (see Section [TV C[) 

<&2 = -~A, R = AA, $ u = ip, (3.3) 

where A = ±(A + + A_) and p = ~(A+ - A_), together 
with $1 = h(E + iB). These electro- vacuum solutions 
are thus of algebraic type D, unless A = which applies 
to a conformally flat Bertotti-Robinson universe and flat 
Minkowski space. Vacuum direct product spacetimes 
(with p = 0) are Minkowski and (anti-)Nariai spaces. For 
the two Plebahski-Hacyan spacetimes, one and only one 
of the 2-spaces is flat. Therefore, $n = ^|A|, so that the 
condition 2<j>n ± 3^2 = is satisfied. 

More details about some of these background space- 
times can be found in the original works [12, EH, E|| , 
reviews 0,01 or 7 e -g-> m [50l [5ll. [55[| . 

B. Type D background spacetimes 

As will be seen in Section IIV CI the gyraton space- 
times are in general of algebraic type II. However, they 



contain a wider subclass of electro-vacuum solutions of 
type D, which can also be naturally regarded as possible 
background geometries. 

Type D electro-vacuum solutions of Einstein's equa- 
tions are known [(| @, EH, H|| H3, H3|- However, their 
forms are usually different from the parametrization of 
the geometry used here. For this reason we will write 
those type D spacetimes, which belong to our subclass 
of the Kundt family, explicitly. All such spacetimes are 
derived in Appendix [UJ here we only summarize the re- 
sults. 

Although these spacetimes have the same curvature 
scalars as in (|3.3I) . they are not, in general, direct-product 
spaces. In particular, they have a lower symmetry than 
the highly symmetric backgrounds discussed above. 

The A+ = case 
(exceptional Plebanski-Hacyan spacetime) 

As shown in Appendix [UJ all type D solutions natu- 
rally split into two cases. For A + = (i.e., A = — p < 0, 
A_ = 2A) we find a generalization of the exceptional 
Plebahski-Hacyan type D electro- vacuum spacetime [43], 
l5ll . |52| . The metric reads 

ds 2 = dx 2 + dy 2 + 2(Ar 2 — L x x — L y y)du 2 
— 2dwdr + 2(0^2; + a y dy)du , 

where Li(u) and ai(u), i = x,y, are arbitrary functions of 
the coordinate u only (i.e., constants on each transverse 
space). This corresponds to the metric (|2.ip with P = 1, 
g — 0, h linear in x, y, and ai independent of x, y. 

For at = it reduces to the exceptional Plebahski- 
Hacyan spacetime. It further reduces to the direct prod- 
uct spacetime (|3.1[) when also both L, vanish. Although 
the functions Li ^ do not enter the curvature scalars 
(|3.3p , the geometry of this spacetime is different from 
that of the direct-product spacetimes (for example, it 
contains another shear-free but non-geodesic null direc- 
tion) . 

Nontrivial coefficients a; can be gauged away using 
the transformation (|2.50|) . However, such a transforma- 
tion generates a nonvanishing metric function g and a 
quadratic dependence of h on x,y. It thus seems that 
the case =/= is indeed a nontrivial generalization of 
the exceptional Plebahski-Hacyan spacetime. 

The A + / case 

In the case when the transverse space has a nonvanish- 
ing curvature A+, the metric of type D electro- vacuum 
solutions is given by the metric functions 

A = |, k = , (3.5) 

h = U+X 2 1 g = 0. (3.6) 
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Here the functions P and Q can be written as 

P=l + - A K+{x 2 +y 2 ) , 

1 4 (3.7) 
Q = q (l- ^A + (x 2 +?/)) +1xx + q y y , 

respectively, where qo(u), q x (u) and q y (u) are constant 
on the transverse space. 

In the case A + > 0, when the transverse space is a 
sphere, the solution for A can also be rewritten as 

A = C [cos 6 cos 0' + sin 6 sin & cos(0 - <t>')} , (3.8) 

where and <fi are standard spherical coordinates, 
cf. (I5.16|) . and C, 8', and <f>' are (possibly u-dependent) 
transverse constants equivalent to qo, q x and q y . 

A slightly more general parametrization of these space- 
times can be found in Appendix [U1 

C. Kundt waves without gyratons (ji — 0) 

Kundt waves on direct product spacetimes 

Now we briefly describe more general Kundt space- 
times of the form (|2.1[) which, however, still do not con- 
tain a gyratonic matter. In such a case, the source func- 
tions ji vanish, i.e., p, q = 0, and the field equation (|2.46|) 
reduces to 

AA + 2A+A = 0. (3.9) 

Let us first consider the trivial solution A = 0; the general 
case is discussed below. 

Since p = 0, it is possible to use the gauge transforma- 
tion (|2.50p to eliminate both k and g, cf. (|2.45[) . Con- 
sequently, we obtain cij = everywhere, and the metric 
simplifies to 

ds 2 = dsl g - 2h(u,x ,y)Au 2 , (3.10) 

where ds 2 g is the metric of direct product spacetimes 
(some of which are listed in Table Q} given by (|3.1[) . 
For nontrivial profile functions h, this class of solutions 
can be interpreted as specific exact Kundt gravitational 
waves which propagate in flat, (anti-)Nariai, Bertotti- 
Robinson, or Plebahski-Hacyan universes (see and, 
for the limit of impulsive waves, HU). 

Indeed, from the corresponding NP scalars (|4.17[) . 
P~2"Tj) (cf. Section IP7C|) . by using (@T7]) it follows that 
(|3.3p remains unchanged, and, in addition, there is 

* 4 = (P 2 h x ) x , $22=P 2 ^, CC -, (3.11) 

where C = (% + M/)/v2- When ^/^ ^ 0, such spacetimes 
are thus of types II or N, and in general, contain a null 
radiation field characterized by $22 • I n particular, pure 
vacuum gravitational waves of this type (which propa- 
gate on a vacuum or electrovacuum background space) 



are given by the condition $22 = 0, so that their profile 
functions h must be of the form 

h = F(u,0+?(u,(), (3.12) 

where T(u, Q) is any function, holomorphic in £. 

Kundt waves on type D backgrounds 

Similarly, we can also describe gravitational waves 
propagating on general type D backgrounds discussed 
above. Indeed, the field equations for A, k, and g are 
linear, and the equation for h is linear in h (with nonlin- 
ear terms with A and k as a 'source'). We can thus easily 
superpose a pure gravitational-wave contribution of the 
form (|3.12p on top of any background metric function h, 
keeping the values of A, K, and g unchanged. In par- 
ticular, considering the exceptional Plebanski-Hacyan 
type D backgroun d (13. 4p . the family of gravitational 
waves described by [591 ] is obtained. 

More general Kundt waves 

The equation (|3.9p is the special Helmholtz equation on 
the transverse space such that the coefficient of the 'mass' 
term is exactly given by the curvature of the transverse 
space. Its general solution can thus be parametrized by 
a single function C(u, Q, holomorphic in £, as 

A = £, c + £ f - 2C (log P) c - 2£ (log P) if . (3.13) 

Again, the functions k and g can be gauged 
away, k = g = 0, and the electro-vacuum condition 
l = p = q = implies A h = 0, see (12.48[) . However, now 
we have an additional contribution to h thanks to a non- 
trivial A, cf. {US]): 

+ (3.14) 

We have thus obtained an explicit form of a general 
Kundt electro- vacuum spacetime (|2.1[) . Apart from non- 
trivial H , these most general gravitational waves within 
our class also have nontrivial metric functions a^, given as 
d£ = —ik + P~ 2 (C + £). We are not aware of a discussion 
of such waves in the literature. 

It should, however, be mentioned that some of 
these solutions have unphysical behavior of the metric 
functions — typical solutions of (|3.13|) and (|3.14p have sin- 
gularities or diverge in transverse directions. They thus 
cannot be interpreted as globally well-behaved gravita- 
tional waves. Nevertheless, some of them can be inter- 
preted as external vacuum solutions around a localized 
matter source, e.g., around a beam of null radiation or 
gyratonic matter. Such solutions will be discussed in 
the next section. Here we only note that they can be 
constructed from given matter sources using the Green 
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functions. They are regular and satisfy vacuum equations 
outside the sources. 

Since the 'mass' term in (|3.9p has a special value, this 
equation also admits globally regular solutions. Regu- 
lar solutions for A are exactly those discussed for the 
type D backgrounds, namely, given by (|3.5p , (|3.7[) (or 
(|3.8p ). The solution for h which leads to the regular ge- 
ometry is given by (|3.14[) with sufficiently smooth T ', e.g., 
when it is quadratic in £. 

D. Gyratons on the flat background 

Our class of solutions also contains, as a subcase, the 
original gyraton on a flat background [l8[. Indeed, for 
a vanishing cosmological constant and electromagnetic 
field absent, the background is Minkowski space. If we 
admit only an r independent gyraton source (i.e., if we 
assume divj = 0, j u = i) and if we employ the gauge 
2 = 0, we immediately obtain the solution discussed in 

IV. PROPERTIES OF THE GYRATON 
SOLUTIONS 

A. Character of the gyratons 

Now we concentrate on nontrivial gyratons contained 
in the above class. A characteristic feature of the gy- 
ratonic matter is a nonvanishing source ji m (|2.5p or, 
equivalently, its two potentials p and q, cf. (|2.43p . The 
gyratonic matter moves with the speed of light, as can 
be identified by inspecting the dependence of the metric 
function on the coordinate u. From the form of the met- 
ric (|2.ip we infer that u is a null coordinate, with null 
generators given by the principal null congruence d r . All 
of the metric functions can depend on this coordinate, 
and this dependence is not restricted by the field equa- 
tions. It means that the profile of the gyraton in the 
u direction can be prescribed arbitrarily. Thanks to a 
trivial r dependence of the fields, such a profile remains 
essentially unaltered (except for the "cooling effect" dis- 
cussed below). This can be understood as a motion of 
the gyraton in the direction of the null congruence d r . 

The characteristic spatial components of the gyraton 
stress-energy tensor represent a possibility of an internal 
energy flow of otherwise null radiation. It can be natu- 
rally split into two components. 

The divergence- free component, controlled by the 
source potential q, corresponds to a "rotational" part of 
the energy flow. However, since the gyratonic matter 
is null, the nature of the "rotation" must be internal — 
it describes a spin of the null fluid. This kind of the 
source was discussed in the context of the gyratons in 
flat spacetime [U and in anti-de Sitter space [26| . 
From the field equation (|2.46p we observe that this "ro- 
tational" part of the source gives rise to the component 



of the metric function dj determined by the potential A 
via (|2.35p . This component is independent of the gauge, 
so the presence of the "rotational" part of the gyraton 
source necessary leads to the nondiagonal component at 
in the metric (|2.ip . 

The rotation-free component of the gyraton source, 
controlled by the potential p, has a different character. 
As can be read out from the conservation law (|2.3ip . 
the source with a nontrivial divergence div j describes an 
internal flow of the energy in the gyraton beam which 
changes its internal energy j u with r. We could thus un- 
derstand the p component of the source as some kind of 
"cooling" which steadily decreases the energy density of 
the gyraton beam. Such a kind of the energy transfer is 
not very plausible physically, mainly because the cooling 
should occur in matter moving with the speed of light. 
It inevitably leads to an unnatural causal behavior of the 
source. 

Indeed, it is easy to check that the gyraton stress- 
energy tensor (|2.5p (composed by either a p or a 
q component) does not satisfy neither a null, weak, 
strong nor dominant energy condition. However, for a 
spinning matter it is not so surprising — bad causal be- 
havior is typical for spinning relativistic objects when 
they are idealized excessively. 

From the equation (|2.45p we also observe that the p 
part of the source controls the combination g — A_k of 
the metric functions. Splitting its influence between g 
and k is just a matter of a gauge choice. We have already 
discussed that it is possible to eliminate either of them 
but not simultaneously. A gyraton source composed just 
from the rotation-free component (p ^ 0, 9 = 0) thus 
does not necessary lead to a nondiagonal component at 
in the metric — its influence can be gauged away entirely 
into the metric function g, and vice versa. 



B. Geometrical properties of the principal null 
congruence 

Let us now briefly discuss geometrical properties of 
the gyraton solutions. The additional property (|2.2p . 
d r ai = 0, characteristic for the subclass of spacetimes dis- 
cussed here, has a consequence that the null vector k is 
recurrent, [g, [23[, namely, 

k a -,p = {-d r H)k a k p . (4.1) 

The null character of k and the condition (|4.ip also im- 
ply that the null congruence with tangent vector k is 
geodesic, expansion- free, sheer-free, and twist-free and 
thus belongs to the Kundt class. 

The condition (|2.2p and the condition (|2.24[) (the func- 
tion H is at most quadratic in r, with a constant coeffi- 
cient in front of r 2 ) guar antee that these Kundt metrics 
are of the CSI type [2{|. For these metrics it was also 
shown that there exists a (u-dependent) diffeomorfhism 
x 1 = x z (u,x k ) such that the transverse metric (|2.9p can 
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be made u-independent. We have already used this prop- 
erty at the very beginning when we applied the conditions 
(|2.12|) . Moreover, the transverse space is locally homo- 
geneous. 

For the complete four-dimensional gyraton spacetimc, 
it was demonstrated in [3(| that there always exists a 
related locally homogeneous spacetime which has invari- 
ants that are identical to those of the Kundt CSI metric. 
This "background" metric can be obtained by setting 
dj = g = h = 0, 



rfs bg = ds^_ — 2 du dr 



A_r 2 du 2 



(4.2) 



which is exactly the metric for direct-product background 
spacetimes (|3.ip . 

The condition 1)2. 2[) is also equivalent to the fact that 
the 2-spaces orthogonal to the transverse spaces are 
surface- forming. 

It is a general property of the Kundt family that the 
vector k is the principal null direction of the spacetime. 
By determining its degeneracy we can thus identify the 
algebraic type. To proceed, it will be convenient to intro- 
duce an aligned complex null tetrad {k, I, to, to}. There 
exists a standard choice of such a tetrad in the context 
of the Kundt family of spacetimes [|| , namely 

k = d r , 

l = 8 u - P 2 (a x d x + a y dy) -(H+ \a 2 ) 8 r , 



P 

to= —(d x + id y ) , 
_ p 

m= -j={d x - id y ) . 



(4.3) 



The spacelike complex vectors m a and m a are tangent to 
the transverse space. Clearly, [k,l] = — (d r H)k, so the 
space spanned on k and I is indeed surface-forming. The 
dual frame in the space of 1-forms reads 

9 (fc) = dr+ (H + ±a 2 )du , 
9 W = du , 

e (i7i) = _L( da; _j d2/ ) + -^={a x ~ia y ) du , ( 4 - 4 ) 
V 2P v 2 

(m) = _^ — (dac+idy) + -^=(a x +ia v ) du . 
V2P 

Calculating the Newman-Penrose spin coefficients 
with respect to this tetrad (see the following section for 
the nontrivial ones), we recover again the general prop- 
erties that the congruence is nonexpanding and non- 
twisting (/3 N p = 0), sheer- free (cr NP = 0), geodesic and 
affinely parametrized (k np = e N p = 0). In addition, from 
k np = 7T NP = £ NP = 0, it follows that the tetrad (|4. 3|) is 
parallelly transported along the null congruence. More- 
over, the condition (|2.2p is directly related to the vanish- 
ing coefficient r NP = 0. 



C. NP formalism in complex coordinates 

It turns out to be more convenient (and common in the 
literature on Kundt spacetimes) to introduce complex co- 
ordinates in the transverse space. Instead of conformally 
flat real coordinates x and y we will now use the complex 
coordinates £ and £ such that 



1 



C = -n= (x + iy) 



(4.5) 



The coordinate 1-forms and vector fields transform as 



dC = -j= {dx + idy) , d c = -^=(d x - idy) 



1 



1 



(4.6) 



d( = —(dx - idy) , = -^=(d x + idy) , 



and the transverse Laplace operator (|2 . 18[) on any scalar 
tp becomes 

Aip — 2P 2 i/> iCC - • (4.7) 

Instead of the real 1-form components a, it is custom- 
ary to introduce a complex function W(u, C, Q by 

W = -a c = -^(a x -ia y ) , W = -a f . (4.8) 

Substituting for the potentials via (|2 . 35[) and using 
= — = — i, = e^q = 0, we find that 

W = -(n + i\) x , W=-(K-iX)^. (4.9) 
We also obtain 

a 2 = 2P 2 WW , 
A K = div a = -P 2 (W x + W x ) , 

A A = -rota = -b = iP 2 f Q c, = iF>2 { W C ■ 

(4.10) 

The metric (|2. 1|) in complex coordinates then reads 

ds 2 = -^dCdC - 2dudr - 2Hdu 2 - 2(Wd(+Wd()du , 

(4.11) 

the canonical form (|2.13p of P is 



P = 1 + -A+CC , 
and the Maxwell tensor (12.41) takes form 



F = Edr A du + B -^d( A d( 

pz 

+ o-£ du A d( + (Jq du A dC 



(4.12) 



(4.13) 



where = (a x — ia y )/^/2. 

The tetrads (|4. 3[) and (|4.4[) are closely related to the 
introduced complex coordinates: 

k = d r , 

l = d u + P 2 (Wd^ + Wd ( ) - (H + P 2 WW)d r , 

TO = P8q, 

to = Pdt^ , 

(4.14) 
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and 



we obtain 



(fe) =dr+(H + P 2 WW)du , 

9 (,) = du , 
6( m ) = i dC - PWAu , 
qM = 1 dC - /Wdu . 

The list of nontrivial NP coefficients is then 

?2t 



(4.15) 



7np = 

a NP = 



\p 2 (w rc + w x ) , 

P(H + P 2 WW) X , 

d r H+\{(P 2 W) x ~(P 2 W) x ) 



-iP- 

2 ■ 



(4.16) 



The source equations can be recovered in the Newman- 
Penrose formalism by comparing the components of the 
Ricci tensor with the corresponding components of the 
electromagnetic and gyraton stress-energy tensor. The 
general form of nonvanishing Ricci scalars for the metric 
(|4.1ip is listed in Appendix IBl in equations (|B2[) . After 
substituting the potentials k, A, cf. equations (|4.9[) . and 
using the relations (|4~8]) . (|4.10|) . we find 



(J)., 



2' 
P 

~2 
1 

2 



~2p« + (.g - A_k) + i(-A A + A+A) 



rA(j- A_k) + iP 2 (a^b x -a^b x ) 
1 6 2 + A + a 2 + gAn + d u An 



R 



+ Ah- 
24A NP 



4A 



(4-17) 

The constants A and p have entered these expressions 
via combinations A = \{A + + A_) and p = ^(A + — A_) 
of the constants A± which parametrize the metric (|4.1ip 
through (12.24j) and (|2.11[) . Their relation to the cosmo- 
logical constant A and the electromagnetic energy den- 
sity p is established by comparing these components to 
the cosmological term and to the corresponding compo- 
nents of the stress-energy tensors. For the electromag- 
netic field, the nonvanishing components are 



(REM 

<T> EM 
v 22 



1 

2 P > 

2P 2 pn X K 



(4.18) 



o, 
p 

1 . 

2 Ju ~ 



(4.19) 



P z 



The first and second source equations (|2.42[1 and (|2.48|) 
are obtained from the above components $12 and $22, 
respectively, by realizing that 



k = (p + «?),c 



3c = (p-k),c ■ 

Finally, in Appendix[B]we also present the Weyl scalars 
(|B3j) in the form which follows directly from the metric 
(|4.11[) without using the field equations. If we introduce 
the potentials and employ the field equations, the non- 
trivial scalars reduce to 



(4.20) 



1 



-A , 



9h = PAk, 



(4.21) 



*4 = r{P 2 (p- iA-\) x ) x +{P 2 h,c) x 

+ [A+A + iA-K + ip + id u ] (P 2 \ x ) c 

+ 2iP 2 q. c ( K + i\). c +p(P 2 K. c ) - 2AP 2 (k. c ) 2 



where h is given by (|2.49|) . 

Since ^0 = ^1 = and, for a nonvanishing cosmolog- 
ical constant, $2 7^ 0, we conclude that the vector k is 
the double degenerate principal null direction and the 
gyraton spacetime is of the algebraical type II. The con- 
ditions for further algebraic degeneracy to the type D 
are, in the vacuum case, discussed in Section IIII Bl and 
in Appendix [Cj In the nonvacuum case these conditions 
are rather strong: for example, there are no nontrivial 
type D gyratons with p = 0. 

For the vanishing cosmological constant A = 0, the 
presence of a nontrivial rotational gyratonic matter 
(given by the potential q; terms with n and p are not 
significant, as they can be cancelled by a suitable gauge) 
guarantees that the spacetime is of type III. The space- 
time reduces to type N only for A = and q = 0. 

Comparing (|Q7|) with P~2l"j) . or directly (|B2|) with 
(|B3|) . we find that $12 and ^3 are closely related, namely 



$12 



(4.22) 



*3 = Pg,t ■ 

The radiative characteristic of the gravitational field ^3 
is thus determined by the matter component $12, up to 
the term which can be controlled by the gauge. In the 
gauge g = 0, we have directly $3 = — $12. 



D. NP components in the gauge invariant tetrad 



with p given by (|2.7[) . Similarly, for the gyratonic matter 



The tetrad ([43]) introduced in Section IIVBI is par- 
allelly transported and transverse-surface forming, how- 
ever it is not gauge invariant under the transformation 
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(|2.50[) . Consequently, the spin coefficients and curva- 
ture scalars can have rather nontrivial dependence on 
the gauge transformation. It is possible to choose an- 
other tetrad which is gauge independent. Actually, it 
is a tetrad which is also well known in the context of 
the Kundt family of solutions Q ■ This is related to the 
tetrad (|4. 3|) by a null rotation with the vector k fixed, 
i.e., it is also aligned with the principal null congruence. 
Explicitly, this null rotation is 



Similarly, for the gyratonic matter we obtain 



k — k , 

m = m + Kk , 
I' = 1 + Km + Km + KKk 



(4.23) 



with the complex parameter K = —PW. Thus, the 
gauge invariant tetrad of null vectors is 



k' = d r , l' = d u - Hd r , 
m' = P{d z - Wd r ) , m' = P(d c - Wd r ) 

and for the dual frame of 1-forms we obtain 



(4.24) 



6 (A:) ' = dr + Wd( + Wd( + Hdu , 9 (i) ' = du , 

(4.25) 

The corresponding nontrivial spin coefficients are 
A^p = , 



ew = -de , e (m) ' = -de 
p s ' p s 



i 



Tnf 



/4 



P(H X - Wd r H - d u W) , 
l -(d r H-Ub) , 



(4.26) 



and the nonvanishing Ricci scalars are 



-2ip\ + (g - A-k) + i(-A A + A+A) 
r A(g- A_k) + 2P 2 (a^g x +a c g^) 



(4.27) 

For the electromagnetic field, the nonvanishing compo- 
nents are 





1 




$11 = 


2 P ' 




P 




$12 = 






~2 






1 




*22 = 


2 


r t 










+ 





-iPp\c i 
^ = 2P 2 p\ c \- c . 



(4.28) 





= 0, 




P 


*!2" 


- J3( 




1 . 


^2 


- 2^ 



(4.29) 



The nonvanishing Weyl scalars read 
P 

tf 3 = -iPAA, c + -(p-ig), c , 

* 4 = K p2 (P - iA - A ),c), c +(^ 2 ^c),c (4-30) 
+ [A_A + iA_k + ip + i<9„] (P 2 A :C ) 
+ 2P 2 p, c ( K + U), c + p(P 2 K C ) + 2A_P 2 (A, C ) 2 . 



E. Electromagnetic field 



In our ansatz made in Section III Al we allowed the 
spacetime to be filled with the electromagnetic field (|2.4[) . 
This field does not have its own dynamical degrees of 
freedom — it is specified just by two constants E and B. 
In the presence of a gyraton, this electromagnetic field 
is modified through the du A aidx 1 terms. However, the 
transverse 1-form cr, is uniquely determined by the gyra- 
ton, see equations (pT34"|) . (pT38)) . 

The Maxwell tensor (|2.4[) can be split into two parts 



F = E(dr Adu + duA dn) + B(P~ 2 dx Ady + duA dX) . 

(4-31) 

It is interesting to observe that the 2-form proportional 
to the constant B is the four-dimensional Hodge dual of 
the 2-form proportional to E. Thus, the Maxwell ten- 
sor has a familiar structure of a linear combination of 
dual 'electric' and 'magnetic' parts. Moreover, after sub- 
stituting (|2.38p for cp, the field equations depend only 
on the 'weights' E and B of the electric and magnetic 
parts through the constant p = (E 2 + B 2 ) (via the 
constants A± = A ± p) . The geometry of the spacetimes 
thus does not depend on a particular splitting of the elec- 
tromagnetic field. 

To inspect the algebraic structure of the electromag- 
netic field, we need the tetrad components $a of the 
Maxwell tensor. With respect to the parallelly trans- 
ported tetrad (|4.14[) . we obtain 



*o = 0, 

$i = ^(E + iB) 



(4.32) 



$ 2 = -P(E + iB)k x , 
while with respect to the gauge invariant tetrad (|4.24|) 
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we get 



*o = 0, 

$1 = \(E + iB) , 
$ 2 =iP(E + iB)X ti 



(4.33) 



It follows that the electromagnetic field is aligned with 
the principal null direction k of the gravitation field, but 
this vector is not a double degenerate vector of the field. 

The corresponding tetrad components of the elec- 
tromagnetic stress-energy tensor = k£ <&a^b have 
been listed in (|4~T5|) and (g^Sl). Notice that the and 
$!" components of (|4. 18[) can be simultaneously can- 
celed by the gauge choice k — 0. This choice also cancels 
the component $2 in (|4.32[) . 



V. GREEN FUNCTIONS 

In Section |TT] we demonstrated that for our ansatz 
the Einstein-Maxwell equations effectively reduce to the 
Poisson equations 



Aip = —s 



(5.1) 



(e.g., equations (|2.41[> . (|2.48|) for A and h), and to the 
Helmholtz-Poisson equations 



A ip + R±ip 



(5.2) 



(equations (|2.46[) and (|2.28|) for A or 6). These equations 
on the two-dimensional transverse space can be solved 
using the Green functions G(o) and G(i), respectively. 
Such functions satisfy 6 

A G(o) (2, x) — — 8{x, x) , (5-3) 
[A + Ri_]G {1) (x, x 1 ) = -S(x, x 1 ) , (5.4) 

where x and x' are points in the transverse space. The 
solutions are then given by the integral over the corre- 
sponding sources 



ij)(x) = J G(„)(x, x') s(x) y g' ± d 2 x' . (5.5) 
In particular, 

X(u,x) — I G(o) {x, x') b{u, x') \J~g^_ d 2 x' , 

(5.6) 



h(u,x)= / G( )(x,x')x(u,x') Jg'±d x' , 



In the case when there exist normalizable zero modes, one has to 
subtract a projector to the space of these modes from the delta 
function on the right-hand side. See a discussion in Appendix [Dl 



with the source 

X = -i-gAA+pAK + 2a ! p )i , (5.7) 

and 

A(u,x) = 2 / G(!)(x, x') q(u, x') J g' ± d 2 x' , 



(5.8) 



b(u,x) = 2 / G(!) (x, x') rot j(u, x') J g' ± d 2 x' 



It follows from the Einstein equations that the two- 
dimensional transverse space is a maximally symmetric 
space of constant curvature R± = 2A + , i.e., a plane E 2 
for A + = 0, a sphere S 2 for A + > 0, and a hyperboloid 
H 2 (Lobachevsky plane) for A + < 0. The corresponding 
Green functions are known explicitly (see, e.g., [60] ) and 
they are discussed in more detail in Appendix |Dj Here 
we present only those results which are important for 
solution of our problem. 



A. Green functions for 2-plane E 2 

For A + = the transverse space is the flat plane. Both 
the Green functions coincide and they have the form 



1 



G< v) (x,x') = -— \og£(x,x') 

Z7T 



(5.9) 



Here £{x, x') is the distance between the points x and x' . 



B. Green functions for 2-hyperboloid H 

For A+ = — 1/L 2 < the transverse space (|2.9|) is hy- 
perboloid of a constant negative curvature, L being the 
curvature radius. It can be parametrized by different use- 
ful coordinate systems. Here, we list some of them which 
are frequently used in the literature, namely, hyperspher- 
ical, Poincare, Lobachevsky, and projective coordinates, 
respectively: 

dsl = L 2 ( dp 2 + sinh 2 p d(f> 2 ) 
L 2 



= -^(d^ + d<) 
= L 2 ( dp 2 + cosh 2 fi dr 2 ) 



(5.10) 



( da; 2 + dy 2 ) 



(l+\A+(x*+y*)y 
Relations of the coordinates to the projective ones are 



x = 2L tanh — cos <fi , y = 2L tanh — sin < 



x = L 



At 



< 2 + (l + zf 



y = L 



2L cosh p sinh r 
cosh p cosh r + 1 



t 2 + z 2 - 1 

t 2 + (l + z) 2 ' 

2L sinh p 
cosh p cosh t + 1 



(5.11) 
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Clearly, the conformally flat coordinates used in the text 
correspond to the projective coordinates with the choice 
(|2.13p and to the Poincare coordinates with the choice 

Because of the maximal symmetry of the space, the 
Green functions can be expressed only in terms of the 
geodesic distance between the points £(x,x') or, more 
conveniently, of its function 

r](x,x') = cosh ^~A+£(x,x')j . (5.12) 

The function r\ in an explicit form reads 

rj = cosh p cosh p — sinh p sinh p' cos(0 — cf> ) 
„, (z-z') 2 
2zz' 

= cosh p cosh p! cosh(r — r ) — sinh p sinh p' 

A + (x 2 +y 2 ) \ f 1 A + (x' 2 +y' 2 ) 



1- 



A + (xx'+yy') 



(l+^M) (l+ 



A + (x' 2 +y' 2 ) 
4 



(5.13) 

Clearly, 77 £ [l,oo), with rj = 1 corresponding to coin- 
cident points, and 77 — ► 00 to an infinite distance. Using 
these quantities, the Green functions in question are 



G {0) (x,x') = -±- log ^-i + l) , 
G (1) (x,x') = -±- ^log L-i + il 2 



When rj — > 00, the Green functions tend to zero. 



C. Green functions for 2-sphere S 2 



(5.14) 



The generic solution for the Green functions on a sphere 
is a linear combination of the Legendre functions Q>(?7) 
and P„(?7). However, the requirement of regularity at 
the antipodal point 77 = —1 singles out their particular 
combination. Also, one has to be cautious since both 
equations (|5.ip and (|5.2p on the compact sphere have 
normalizable zero modes — see Appendix [Dl for more de- 
tails. Eventually, the Green functions read 



G(q)(x,x') 



1 

47T 



log (1 - rj) 



G {1) {x,x') = -^(r 1 \og{l-rj) + l) 



(5.19) 



Moreover, the left-hand side of (|2.48p is the Lapla- 
cian defined on a compact sphere. The integral of the 
Laplacian over the sphere has to be zero. This property 
imposes an integral condition on physically acceptable 
distributions of the stress-energy tensor, namely 



/ xVgId 2x = 

Js 2 



(5.20) 



Similarly, for the equation (|2.28|) we also get an integral 
constraint 



/ cos 9 rot j y/g± d 2 x 
Js 2 







(5.21) 



Because of this property the zero modes do not contribute 
to the components of the metric. The constraints 
(|5.20[) and (|5.2ip appear only because S 2 is compact and 
are analogous to the property that closed worlds must 
have zero total energy, charge, or angular momentum 

MM. 



When A + = l/L 2 > the metric ([23]) describes a 
sphere 



ds\ = L 2 (d6» 2 +sin 2 ^ 2 



dx 2 + dy 2 



(l+lA + (x 2 +y 2 )Y 

(5.15) 

Spherical coordinates (9, (f>) are related to the projective 
coordinates (x,y), used in (|2.ip with P given by (|2.13p . 
via the coordinate transformation 

9 9 
x = 2L tan — cos <f> , y = 2£ tan — sin <fr ■ (5.16) 

Then the Green functions are functions of 

r](x, x') = cos (y/X^e(x, x')) . (5.17) 
Here r\ varies in the interval [—1, 1] and has the form 
rj = cos 9 cos & + sin 9 sin 9' cos(</> — (/)') 



+ A+{xx'+yy') 



^ | A + (x 2 +y 2 ) ^ ^ | A+{x' 2 +y' 2 ) \ 



VI. CONCLUSION 

We presented a new class of gyraton solutions on 
electro- vacuum background spacetimes which are formed 
by a direct product of two constant-curvature 2-spaces. 
These involve the (anti-)Nariai, Bertotti-Robinson, and 
Plebahski-Hacyan spacetimes in four dimensions. The 
background geometries are solutions of the Einstein- 
Maxwell equations corresponding to the uniform back- 
ground electric and magnetic fields. The gyraton solu- 
tions are of Petrov type II and belong to the Kundt fam- 
ily of shear- free and twist-free nonexpanding spacetimes. 

Gyratons describe the gravitational field created by a 
stress-energy tensor of a spinning (circularly polarized) 
high-frequency beam of electromagnetic radiation, neu- 
trino, or any other massless fields. They also provide 
a good approximation for the gravitational field of a 
beam of ultrarelativistic particles with a spin. The gy- 
ratons generalize standard pp -waves or Kundt waves by 
admitting a nonzero angular momentum of the source. 
This leads to other nontrivial components of the Ein- 

xT U i, in addition 



(5.18) stein equations, namely G u i + Ag n 
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to the pure radiation wit-component which appears for 
pp -waves or Kundt waves. 

We have shown that all of the Einstein-Maxwell equa- 
tions can be solved exactly for any distribution of the 
matter sources (see Section [II Fl for a summary), and the 
problem has been reduced to finding the scalar Green 
functions on a two-dimensional sphere, plane or hyper- 
boloid. These Green functions have been presented in 
detail in Section |Vl Special cases of these gyraton so- 
lutions and their properties are discussed in Sections IIIII 
and [IVl 

We have also studied the gyraton solutions using the 
Newman-Penrose formalism. The characteristic term ai, 
describing the rotational part of the gyraton, generates 
the nontrivial Ricci $12 and Weyl ^3 scalars, in addition 
to the case of pure pp and Kundt waves. Curiously, there 
exists a very simple relation (|4.22p between them. 

To complete our investigation, we have also studied 
gyratons on more general type-D backgrounds (includ- 
ing the exceptional Plebahski-Hacyan spacetime) which 
are not direct-product spaces. In addition, in Section 
MI CI we have identified a special subclass of the gyra- 
ton solutions — general vacuum Kundt waves which also 
contain cases previously not discussed in the literature. 

A natural next step would be the study of gyratons 
in a full family of Kundt spacetimes, especially on con- 
formally flat backgrounds, including the (anti-)de Sitter 
universe. Another generalization could be their exten- 
sion to higher dimensions, where, however, one has to 
deal with a richer possible structure of the transverse ge- 
ometries. 
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Appendix A: THE EINSTEIN EQUATIONS 



The inverse to the metric (12. ip is 
grdpdu = P 2 {OA + d y dy) - 2d u d r (Al) 
+ 2P 2 (a x d x + a y P 2 dy)d r + 2(H + ^a 2 ) d r d r . 

The stress-energy tensor of the electromagnetic field 
T EM , according to the definition 

= e (f^F vX g KX - \ 9fiv F kX F kX ) , (A2) 
has nonzero components 



— 2Hp + xe ((T — Ea) 2 , 
xT™ = xe {^{E 2 ~B 2 )a x ~ EB a y - Ea x + Ba y 

xT™ = xe (^(E 2 ~B 2 ) a y + EBa x - Ea y - Ba 3 

rpEM P 

ni xx — p2 ' 



(A3) 



TiEM P 

yy ~ p2 ' 



where the density p = ^^(E 2 +B 2 ) was defined in (|2.7[) . 
The Einstein tensor for the metric (|2.1[) reads 

G ur — A log P , 

G uu = \b 2 + AH+ {d 2 H)a 2 + 2a t d r H. l 



G,, 



1. 



(d r H) div a + <9„div a + 2HA log P 
-b.y - a x (A logP - d 2 H) + d r H, x , 



1 



(A4) 



b, x - a v (A logP - d 2 H) + d r H iV , 



G xx — 9 2 H , 



r — — f) 2 tj 

^yy — p2 »• ' 



where 



A log P — P(P XX + P tVV ) ~ (P% + P 2 y ) , (A5) 



cf. I|2.10p . Here we have used only the metric (|2.ip . with- 
out assuming any other information about the metric 
functions. In particular, we have not used the field equa- 
tions. To be more precise, in the components G U i and 
Gij we employed the fact that P is u-independent; cf. 
relation (|2.12p . However, as we already mentioned in 
section III B[ such a choice is always possible provided 
that the transverse spaces have the same homogeneous 
geometry — which can be derived just from the compo- 
nent G ur . 



Here we present geometric quantities which appear in 
the Einstein field equations, namely the Einstein tensor 
of the general metric (|2.ip and the electromagnetic stress- 
energy tensor corresponding to the field 



Appendix B: NEWMAN-PENROSE QUANTITIES 

In this appendix we present the Ricci and Weyl curva- 
ture scalars corresponding to our subclass of the Kundt 
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family given by the metric (|4.1ip . 



The gauge invariant tetrad 



ds 2 = -^d(dC - 2dudr - 2Hdu 2 - 2(W&C+WdC)&u , 

(Bl) 

with functions W(u, £, C) an d -P(CjC) independent of r, 
or r and u, respectively. The expressions below have not 
employed any field equations. 



The standard tetrad 



With respect to the tetrad (|4.14[) the only nonvanishing 
Ricci scalars are 



$11 
$12 

$22 



1 
4 
P 
~2 

p2 

2 



(d 2 r H + A\ogP) , 
d r H x + (A logP)W - -ib t f 

J?- + 2H < z-i(Wbz-Wb, ( ) 



+ 2(MogP)WW - (d r H+d u )(W x +W x ) 
R = 24A NP = 2(-d 2 H + A logP) , 
and nonvanishing Weyl scalars read 

1 



(B2) 



y 2 = -(d 2 H-AlogP), 
o 



*3 = 



P 



d r H x - (AlogP)VK- -ib x 



(B3) 



* 4 = -iP z Wb x - (MogP)P 2 W z - (P'd u W) x 
+ (P 2 H X ) X ~ (d r H)(P 2 W) x . 



Here, A logP is given in (|A~5j) , and b = iP 2 (W x - W x ) 
as in (|4.10p . The components of the electromagnetic field 
are 



$n = 0, 



$i = -(E + iB), 



(B4) 



$2 = P{-(?i +iBW) 
where cr^ = -i=(cr x - ia y ). 



With respect to the gauge invariant tetrad (|4.24[) we 
obtain that the only nonvanishing Ricci scalars are 



$11 = 


\(d 2 r H + AlogP) , 




$12 = 


P 
~2 


-W{d 2 r H) + d r H x - -ib x 


} 


$22 = 


P 2 
2 


\ b 2 — 

— + 2P cc - - 2(Wd r H x +Wd r H x ) 




+ 2{d 2 H)WW - {d r H+ d u )(W x +W x ) , 



(B5) 



the nonvanishing Weyl scalars read 

* 2 = i(d r 2 P-AlogP), 
6 



$3 = 



p 



~W{d 2 H) + d r H x - -ib x 



(B6) 



* 4 = P 2 W 2 (d 2 H) - {P 2 d u W) x + (P 2 H X ) X 
- P 2 W(d r H) x - {P 2 Wd r H) x , 

and, finally, the scalars for electromagnetic field are 

$o = 0, 
1 



$i = -(E + iB) , 



(B7) 



$ 2 = -P(EW + a ( ) . 

Appendix C: ALL SPACETIMES OF TYPE D 

The purpose of this appendix is to derive all electro- 
vacuum solutions of the algebraic type D within the class 
considered. Inspecting the field equations (|2.8ip - (|2.83|) 
without the gyratonic matter, i — p = q = 0, we find 



9 = A_k , 
P 2 A c£ + A+A = , 



\CC 

p 2 h 



cc 



o, 



(CI) 
(C2) 
(C3) 



with h given by (|2.84|) . Using the gauge transformation 
(|2.50p . we could eliminate both k and g, but this is not 
necessary in the following. 

From (|4.2ip we infer that the vector k = d r is a double 
degenerate principal null direction. When ^ 7^ 0, i.e., 
for a nonvanishing cosmological constant, the condition 
that there exists another degenerate null direction is 



cf. 



3$4$ 2 

This reduces to the relation 



= irA-(P z \. 



{P 2 h x ) c + [A+A + A_k + id u ] {P 2 \ C ) X 



(C4) 



(C5) 
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Taking into account the r dependence, we obtain the fol- 
lowing two conditions : 



(^,c), c = o, 

{P% c ) = . 



(C6) 
(C7) 



These equations must be accompanied by the condition 
(|2.1ip for the metric function P, namely, 



P 2 (logP 2 ), c - = A 



The case A+ / 



(C8) 



Integrating (|C2[) . we obtain 
A,< = 4 



p2 . 



(C9) 



where the arbitrary function C(u, () can depend only on 
t and u. First assuming A + 7^ 0, we can substitute the Q 
derivative of (IC9D into (IC2D. which leads to 



A 



£ f -P 2 



A+P 2 



(CIO) 



Now we have to check the consistency of this general so- 
lution for A with (|(J9[) (since we have used the equation 
(|C2[) to obtain A) . It turns out that the solution is consis- 
tent only thanks to the fact that the constants A + which 
appears in (|C2[) and (|C8|) are the same. 

However, the function C in (|C10p is not arbitrary. The 
last condition which must be satisfied is that A is real. 
It is not straightforward to find the consequences for C 
explicitly in a general case. But we can use the freedom 
in transverse diffeomorfhism to transform the solution of 
(jC8[) into a particular form. It will be useful to use such 
transverse coordinates for which P is linear both in C, 
and 



Explicitly, 



P = Pa + Pit + Pit + P2tt 



(Cll) 



(C12) 



with po and pi real constants, and pi £ C, satisfying (as 
a consequence of (|CSJ>) the relation 



P0P2 - P1P1 = ^ A + • 



The solutions (|2. 13[) and (|2.14j) are particular examples 
of such a choice. 

Assuming (|Clip . we can now easily find the reality 
condition for A given by (jClOj) . Since P is real, it requires 



that Q = PA = (2£P C - - C^P)/A + is also real. Taking 
the derivative Qxc> we nn( ^ that it vanishes, and Q is 
thus linear in both £ and 



Q = Qo + Qit + 91C + 92CC • 



(C14) 



Here go and qi must be real and q\ complex transverse 
constants, but they can by it-dependent. The field equa- 
tion (|C2[) now gives the restriction 



Po92 + P2qo =piqi+ piii ■ 

Substituting in (|C9|) , we find 

£ = (Poll - Q0P1) + (PiQ2 - P2qi)t 2 
+ {PoQ2 -P2qo-PiQi +Piqi)t ■ 

In particular, for P given by (|2.13|) we get 



(C15) 



(C16) 



(C17) 



?0 (l--A+CC)+9iC + ?iC, (C18) 



which is equivalent to (|3.7p in real coordinates x and 
y. For P = \J — A_|_ x, which is a generic example of the 
choice (|2.14p . we have Q = Cyf— A+ y and 



(C19) 



, . t-t y 

A = — ic = = c — 

C + C x 



Finally, we should solve the equations (|C3[) and (|C7[) 
for h. As for A, we find that h x = P~ 2 7Y, with H x = 0. 
By substituting into (|C3p . this leads to 



= P 2 h ^ = H x - H(logP 2 



(C20) 



Taking the derivative with respect to t and using (|C8p . 
we find 



(C21) 



= W(logP 2 ), Cf -=^ 



For A_|_ / Owe thus obtain that the conditions (|C3[) and 
(IC7P admit only the trivial solution H — 0. Thanks to 
(I2~i9"l) , for ft we obtain 



h = ^A+A 2 - ^A_k 2 - d u n 



(C22) 



(CI 3) in which k can be set to zero by a proper gauge. 



The case Aj 







If the transverse space is flat, i.e. A + = 0, we naturally 
chose Cartesian transverse coordinates for which P = 1. 
We thus immediately get the conditions 



The case A_ = can be easily discussed separately. 



\(C = 



A 







(C23) 
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which imply A,£ = q\ and 



2-hyperboloid H 2 



A = qo + qiC + 9i C = <7o + fez + O^y 



(C24) 



where Qo(u),q x (u) and are real and qi(u) complex 

functions of u only. The metric 1-form Oj is then 



(k + q y x - q x y),i 



(C25) 



Thanks to such a very special form (|C24[) , it is a gradient, 
and therefore it could be transformed away by a suitable 
choice of k using the gauge transformation. However, it 
would generate a nontrivial function g and contributions 
to the equation for h. Therefore, such a gauge fixing may 
not be the best choice. 

Since in this case h satisfies the same equation as A, 
we can write 

h = Lq + Li( + Li( = Lq + L x x + L y y , (C26) 

with Lq(u), L x (u) and L y (u) real and L\{u) complex 
transverse constants. With the gauge k = g = 0, the re- 
lation (pm?)) gives 



h = L, 



L y y 



(C27) 



where the constant Lq has been eliminated by rescal- 
ing the coordinate u, i.e., incorporating the gauge trans- 
formation (|2.70p . We have thus obtained the metric 
p.4p which is the a* ^ generalization of the exceptional 
Plebahski-Hacyan spacetime. 



Appendix D: SCALAR GREEN FUNCTIONS ON 
CONSTANT CURVATURE SPACES 



The Green functions for the Poisson and Helmholtz- 
Poisson equations (|5.ip and (|5.2p discussed in Section [V] 
are special cases of the Green function satisfying 



[A -£R±\ G(x,x') = -5{x,x') 



(Dl) 



with £ = and £ = — 1, respectively. In our case, the 
scalar curvature R± = 2A + is constant, and the mani- 
fold is either a 2-plane, a 2-sphere or a 2-hyperboloid, 
depending on the sign of A + . Because of the symme- 
try of these spaces the Green function G(x,x') is only 
the function of the geodesic distance £(x, x') between the 
points x and x' . In this appendix we present the Green 
functions for an arbitrary parameter £. 

Since the case of flat plane is trivial, we will discuss 
only the cases of a hyperboloid and a sphere. Let us note 
that on a sphere, for particular values £ = — ^1(1 + 1) 
with integer I > the equation (|Dip has to be modified 
in order to take into account normalizable zero modes. 
Namely, it is neccessary to subtract from the 5-function 
a projector to the space of zero modes. 



For A + < 0, the Green functions can be ex pre ssed in 
terms of the Legendre functions Q v (see, e.g., |60 |) 



G [v) (x,x') = — Q v {rj) , 



where 



-2£ 



(D2) 



(D3) 



and r](x,x') is given by the equations (|5.12p . (|-5. L3[) . 
When the separation between the two points increases, 
?y — > oo and the Green function (|D2|) tends to zero. This 
is a proper boundary condition for Green functions on 
hyperbolic spaces. In our paper only the v = 0, 1 cases 
are important, and the corresponding Green functions 
can be expressed in terms of elementary functions: 



G {0) (x,x') 



1 

2^ 



G(i){x,x) = — Qi{rf) 



47T V 7y 



4^(" l0g (' 



1 



(D4) 



(D5) 



2-sphere S 2 

When A+ > 0, the Green function (|D1| for arbitrary £ 
can be expressed in terms of the Legendre function P„: 



G( v ){x,x') 



1 



4sin7r^ 



(D6) 



see [60|, with v given by (|D3|) and r)(x, x') by the equa- 
tions (|5.17p and (|5.18p . The Legendre functions P^, 
are defined on the interval [—1,1] and are related to 
the Legendre functions P v , Q v on the complex plane as 
F> (77) = ±[P v (r]+i0) +P v (rj-i0)], and similarly for Qs 
[63(. When < 77 < 1 we can observe from the second 
line of (|D6p that there is the required singularity at 77 = 1, 
when both points coincide. For —1 < 77 < 0, the first line 
gives the regularity of the Green function at 77 = — 1 cor- 
responding to antipodal points. 

If v is an integer number, one has to be cautious be- 
cause of a contribution of normalizable zero modes. The 
zero modes must be removed from the spectrum. This 
case thus must be treated separately. 

When v = 0, the corresponding Green function then 
satisfies the equation 



A_ 



AG i()) (x,x , )^~d(x,x') + ^ 



(D7) 



where = f s2 ^/g±d 2 x is the total area of the compact 
transverse space (|2.9p . The term subtracted from the 
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delta function is exactly the projector on the zero mode, 
which for v = is constant. 

The solution of this problem, which is regular at an- 
tipodal points (rj = —1) reads 

G i0) {x,x') = -^-\og(l-r ) ) , (D8) 

with r\ given again by (|5.17p . (|5.18p . One can shift this 
Green function by a constant to satisfy different bound- 
ary conditions at the antipodal points, but this constant 
is unimportant since it will drop out of the equation 
(|2.33p anyway due to the constraint (I5.20|) . 

In the case of the equation (|2.28|) we have v = 1. In 
this case the projector on zero modes is proportional to 
Px(r)(x, x')) = T](x,x') and the equation for the Green 
function reads 

3A 

[A + 2A+] G (i) (x, x') = -S(x, x') + — — rj(x, x') . 

w 471" 

(D9) 

The solution can be found explicitly as 

G (1) (ar,ar') = -^(7 ? log(l-r 7 ) + l) . (D10) 

Similarly to the v = case, one can add to the Green 
function the projector to the zero modes with an arbi- 



trary coefficient to match different boundary conditions 
at the antipodal point. But the component T^f 1 — = 
£i fc <Z,fc of the stress-energy tensor of the gyraton has to 
satisfy the integral constraint similar to (|5.20[) . namely 

/ cos6> q^gld 2 x = . (Dll) 
Js 2 

Because of this property, zero modes do not contribute 
to the components <ij of the metric. Constraints (|5.20[) . 
(|5.21[) appear only because S 2 is a compact manifold. 

Let us finally note that in the case of a general integer 
v = I the modified Green function equation is 

[A+/(/+l)A+] G (l) (x,x') 

21+1 (D12) 
= -6(x,x')+A + ^P l ( V (x,x')). 

47T 

Its solution can be written as 

G (l) (x ) x') = -±-(p l (r 1 )log(l-r,) + U l (r,)) , (D13) 

with Ui (77) being a specific polynomial in rj of the (I— l)-th 
order. In particular, Uq = 0, U± = 1, U2 = \f] + |, 
U3 = |?7 2 + §»?-§, etc. 
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